home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dgecon.z / dgecon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          DDDDGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DGECON - estimate the reciprocal of the condition number of a general
  10.      real matrix A, in either the 1-norm or the infinity-norm, using the LU
  11.      factorization computed by DGETRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, IWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          INTEGER        IWORK( * )
  23.  
  24.          DOUBLE         PRECISION A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      DGECON estimates the reciprocal of the condition number of a general real
  28.      matrix A, in either the 1-norm or the infinity-norm, using the LU
  29.      factorization computed by DGETRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as
  33.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      NORM    (input) CHARACTER*1
  38.              Specifies whether the 1-norm condition number or the infinity-
  39.              norm condition number is required:
  40.              = '1' or 'O':  1-norm;
  41.              = 'I':         Infinity-norm.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) DOUBLE PRECISION array, dimension (LDA,N)
  47.              The factors L and U from the factorization A = P*L*U as computed
  48.              by DGETRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      ANORM   (input) DOUBLE PRECISION
  54.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  55.              NORM = 'I', the infinity-norm of the original matrix A.
  56.  
  57.      RCOND   (output) DOUBLE PRECISION
  58.              The reciprocal of the condition number of the matrix A, computed
  59.              as RCOND = 1/(norm(A) * norm(inv(A))).
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          DDDDGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace) DOUBLE PRECISION array, dimension (4*N)
  75.  
  76.      IWORK   (workspace) INTEGER array, dimension (N)
  77.  
  78.      INFO    (output) INTEGER
  79.              = 0:  successful exit
  80.              < 0:  if INFO = -i, the i-th argument had an illegal value
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.